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Abstract 

We present a technique for derandomising large deviation bounds of functions on the unitary 
group. We replace the Haar distribution with a pseudo-random distribution, a fc-design. Re- 
designs have the first k moments equal to those of the Haar distribution. The advantage of this 
is that (approximate) fc-designs can be implemented efficiently, whereas Haar random unitaries 
cannot. We find large deviation bounds for unitaries chosen from a fc-design and then illustrate 
this general technique with three applications. We first show that the von Neumann entropy of 
a pseudo-random state is almost maximal. Then we show that, if the dynamics of the universe 
produces a /c-design, then suitably sized subsystems will be in the canonical state, as predicted by 
statistical mechanics. Finally we show that pseudo-random states are useless for measurement 
based quantum computation. 

1 Introduction 

There are many results in quantum information theory that show generic properties of states or 
unitaries (e.g. [El IS]). Often, these results say that, with high probabihty, a random state or 
unitary has some property, for example high entropy. However, simple parameter counting shows 
that random unitaries cannot be obtained efficiently. This limits the usefulness of such results since 
no physical systems will behave truly randomly. To make such results more physically relevant, 
it would be desirable to show that these properties are generic properties of unitaries from some 
natural distribution that can be implemented efficiently. Only then could we conclude that we 
would expect to see such properties in natural systems. 

In many cases, the generic properties of unitaries are desirable but randomised constructions given 
by the large deviation bounds are inefficient. We would like to come up with distributions which 
can be implemented efficiently that have similar generic properties. Two examples where the 
best known constructions are inefficient randomised ones are oo-norm randomising maps [15^ [3] 
and locking of classical correlations [9^, fT5]. The results show that, with some non-zero probability, 
random unitaries have the required property. However, there are no known efficient constructions of 
unitaries with these properties. If, on the other hand, we could show that unitaries drawn randomly 
from a set that can be implemented efficiently have the property with non-zero probability, we could 
move an important step closer to finding efficient constructions. (It would not actually provide an 
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efficient construction unless we could find an efficient sampling method.) In fact, for the case of 
oo-norm randomisation, this was done by Aubrun in [3]. 

By random unitaries, we mean unitary matrices distributed according to the unitarily invariant 
Haar measure. In this paper, we will consider replacing the Haar measure with a /c-design. A k- 
design is an ensemble of unitaries such that the k^^ moments are the same as for the Haar measure 
[8]. (A:-designs are formally defined in Section [2]). In particular, this means that the expectation of 
a polynomial in the elements of the unitary matrices of degree at most k is the same whether the 
distribution is the Haar measure or a /c-design. 

The reason for using A:-designs is two-fold. Firstly, because the first k moments are the same we 
would expect similar (although weaker) measure concentration results. Secondly, for k = poly(n) 
(when the design is on n qubits), we might expect to be able to implement the fc-design efficiently 
(i.e. in poly(n) time). Indeed, for k = 0(n/ log n), [13] provides an efficient fc-design constructioru. 

Not only can A;-designs be constructed efficiently, they may even be the product of generic dynamics. 
In [TTj , it is shown that random quantum circuits quickly converge to a 2-design for a quite general 
model of such circuits. It is also conjectured in [14] that random circuits give fc-designs for k > 2 
and k = poly(n) in polynomial time. If a physical system can be accurately modelled by a random 
circuit then, assuming this conjecture, the naturally occurring states will be fc-designs rather than 
fully random states. 

We now summarise some related results in this area. Smith and Leung [22j and Dahlsten and 
Plenio [7j found large deviation bounds for stabiliser states. They showed that, in certain regimes, 
stabiliser states are very likely to have large entanglement. Stabiliser states are 2-designs so our 
results can be seen as a generalisation of this to fc-designs for k > 2 and to other problems. There 
are also some recent classical results related to the present work. Alon and Nussboim [Ij consider 
replacing full randomness with /c-wise independence, a classical analogue of fc-designs, in random 
graph theory. They show that A:-wise independent random graphs with k = poly(log N) {N is the 
number of vertices) have similar generic properties to fully random graphs. 

1.1 Introductory Problem: Entanglement of a 2-design 

We now illustrate our main idea by showing a large deviation bound for the entanglement of a 
2-design, but in a different way to [22\ [7]. 

It has been known for a long time that random states are highly entangled across any bipartition 
|18lll0t [20j. Further, in |16j . it is shown that random unitaries generate almost maximally entangled 
states with high probability. However, generating random states is inefficient so it is an interesting 
question to ask if random efficiently obtainable states are highly entangled. 

Let the system he Ti. = Tis ®'He-, where we label the two systems S and E. Let the dimensions 
be ds and (Ie and d = dsds- Let the overall initial state be any fixed state /jq. Then consider 
applying a random unitary U to SE to get the state •i/' = UpoU^ . Then the von Neumann entropy 
S{'ijjs) of the reduced state tps = trE'4' is close to log2 c?5 (the maximal) with high probability: 

Theorem 1.1 ([16] Theorem 3.3). Let dE ^ ds > 3. Then for unitaries chosen from the Haar 

^Only when we allow for approximate designs, however, we can make the approximation good enough to not 
significantly affect the results. 
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measure 

P(5(V5) < log2 ds-a-(3)< exp ( JAz^)^ ] (1.1) 



(log2 ds)'^ 



whereC=^ and (3 = ^^. 



Now, consider choosing the unitary from a 2-design instead. Later on (Lemma 14. ip . we show that 
Etr-f/^l = "^^^^^ =: /i. Since purity is a polynomial of degree 2, it does not matter if we take the 
expectation over the Haar measure or the 2-design. We now apply Markov's inequality: 



(trV'|>^7) < 



EtrVl 
1 



7 

Using the bound S{'ips) ^ ~log2tr^/'| and some manipulations (the details are in Section this 
can be written as 

P(S(^5)<log2t^s-a-/?)<2-" (1.2) 

where /3 is as in Theorem 11.11 This bound is much weaker than the bound in Theorem 11.11 and, 
in particular, does not show strong concentration as d increases. Later in the paper, we will show 
that choosing unitaries from a /c-design with larger k will give a much stronger bound that does 
give sharp concentration results for large d. 



1.2 Main Results 



We will now state our main results. In the remainder of the paper we will use the following notation 
to identify the distribution we are using. When the unitaries are chosen from a distribution i/, we 
will write to mean Fu^^, the probability when U is chosen from u. Similarly for Ej^, the 
expectation. Usually z/ will be a fc-design. When the distribution is the Haar measure, we will 
write a subscript H. So for the Haar average we write E// for E^/^^^^). 

Our most general result is: 

Theorem 1.2. Let f be a polynomial of degree K. Let f{U) = ^•ajMj(C/) where Mi(U) are 
monomials and let a{f) = \oii\. Suppose that f has probability concentration 

PH(|/-/i| >(^) <Ce-"'5' (1.3) 
and let v be an e-approximate unitary k-design. Then 

ip^d/ - ^ ^) ^ ^ (7)"^ + i + (1-4) 

for integer m with 2mK < k. 

We therefore take a bound for Haar random unitaries of the form Eqn. ll.3] and turn it into a bound 
for fc-designs. For our definition of e-approximate designs, see Section [2j Often, we will use Levy's 
Lemma (Lemma 13. 2p to give the initial concentration bound in Eqn. II. 3[ In this case, a = Q{d) 
(provided the Lipschitz constant (see later) is constant). 

We then apply this to entropy, as a generalisation of Section 11.11 We go via the 2-norm since the 
entropy function is not a polynomial. We find 
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Theorem 1.3. Let u be a 4~" -approximate unitary j^rj^ design on dimension 2" with n > 19. 



Let d = dsdE and 2 < ds < 2"/^° and a > 2. Then 



\{S{i;s)<log2ds-a-P)<8exp2{-^ (? + «) 1 (l-^) 

80 log2 n V 5 



where (3 = j^^. 



We choose a /c-design for A; = -i^q^"^ ^ since this is (up to constants) the largest k for which we have 
an efficient unitary fc-design construction (see Section 12. ip . 

We then move on to apply our results to ideas in statistical mechanics from [T9|. In this paper, 
the authors show that, for almost all pure states of the universe, any subsystem is very close to 
the canonical state obtained from the principle of equal a priori probabilities. This requires the 
dynamics of the universe to produce a random unitary which would take exponential time in the 
size of the universe. We show that the random unitary can be replaced by a A:-design: 

Theorem 1.4. Let Qs the canonical state of the system and ps be the state after choosing a 
unitary from an e-approximate k-design. Let dji be the dimension of the universe's Hilbert space 
subject to the arbitrary constraint R (normally this will be a total energy constraint). Then for 

e<um"\k<'^ 



^u{\\ps-^s\\i>5)<q[j^^ . (1.6) 

Finally, we use results from |llj to show that most states in an 0(l)-approximate state n^-design on 
n qubits are useless for measurement based quantum computing, in the sense that any computation 
using such states could be simulated efficiently on a classical computer. We do this, following [llj, 
by showing that the states are so entangled that the measurement outcomes are essentially random. 

1.3 Optimality of Results 

An important question is how close our results are to optimal, in terms of their scaling with 
dimension d. In Theorem 11.21 we will normally have a = Q{d) so for m constant, we obtain 
polynomial bounds, rather than the exponential bounds for full randomness. This is to be expected: 

Theorem 1.5. Let u be an e-approximate unitary k-design. Suppose also that it is discrete i.e. con- 
tains a finite number S of unitaries. Let f{U) be any function on matrix elements of U and fi be 
any constant. Then either f{U) = fi for all U in ly or for some 5 > 

F^(|/-/i| > 5) >p™„ (1.7) 

where Pmin is the probability of choosing the least probable unitary from v. If the probability is 
uniform, pmin = '^/S. 

Proof. There exists at least one U such that \f{U) — p\ > 5 for some 5 > 0; the probability of 
selecting one is at least Pmin- 
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Corollary 1.6. Our results are polynomially related to the optimal (i.e. the optimal bounds can be 
obtained by raising ours to a constant power). 



Proof. Our results apply for any design, so must obey the bound in Theorem 11.51 for all designs. 
The unitary design construction we use (Lemma 12. 7p has Pmin = d~^^^^ hence the bounds cannot 
scale better than this. □ 

We can also almost recover the tail bound for full randomness in Theorem II. 2[ Suppose for sim- 
plicity that we have an exact design (i.e. e = 0), so that 

( 777 \ m 
^) • 

The optimal m is a5'^/e, which gives 

F^(|/-/i| >S)< Ce-""^'/'. 

So our result allows us to interpolate from Markov's inequality, which gives weak bounds, all the 
way to full Haar randomness and is within a polynomial correction of optimal for the full range. 

The remainder of the paper is organised as follows. In Section [2] we formally define /c-designs 
and what we mean by approximate designs. Then in Section [3] we present our main technique 
for finding large deviation bounds for /c-designs. We then apply this to entropy in Section HI to 
ideas in statistical mechanics in Section [5] and to using /c-designs for measurement based quantum 
computing in Section [H We then conclude in Section [71 

2 /c-designs 

Here we formally define /c-designs. 

Definition 2.1. Let u be a distribution on the unitary group, u is a unitary k-design if 



E 



H 



(2.1) 



for all d^ X d^ complex matrices p (not necessarily valid states). 



We can write this as an equivalent, and for our purposes more useful, definition in terms of mono- 
mials of the elements of the matrices. We will first define what we mean by degree of a monomial 
(or polynomial): 

Definition 2.2. A monomial in elements of a matrix U is of degree {ki,k2) if it contains ki 
conjugated elements and /c2 unconjugated elements. We call it balanced if ki = k2 and will simply 
say a balanced monomial has degree k if it is degree {k,k). A polynomial is of degree k if it is a 
sum of balanced monomials of degree at most k. 

So that, in this definition, UijU*^ is a balanced monomial of degree (1, 1) and UijUki is a monomial 
of degree (2, 0) and is unbalanced. We now state an equivalent definition of unitary A;-designs in 
terms of monomials: 
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Definition 2.3. v is a unitary k-design if, for all balanced monomials M of degree k, 

E^M{U) =EhM{U). 



(2.2) 



That definitions 12.11 and 12.31 are equivalent can be seen by considering matrices p of the form 
\ii,i2, ■ ■ ■ ,ik){ji^j2, ■ ■ ■ ,jk\ in Definition 12.11 Then each element of U®^p{U'^^'^^ is a balanced 
monomial of degree k. Further, each balanced monomial appears for some choice of p. 

We will use state /c-designs, which are related to unitary fc-designs although less general: 

Definition 2.4. Let u he a distribution on states and let uh he the uniform distribution on states, 
which can he thought of as a random unitary being applied to any fixed state. Then u is a state 
k-design if 



E 



(IV')(V'I)' 



E 



mm' 



(2.3) 



By considering unitaries acting on a fixed state, it can be seen that a unitary /c-design can provide 
a state fc-design, although the reverse is not necessarily true. 



2.1 Approximate A;-designs 



There are no known efficient constructions of exact unitary fc-designs. However, for our purposes, 
only approximate designs are required. In Ambainis and Emerson define an e-approximate state 
/c-design using the oo-norm: 



Definition 2.5 ([2j). u is an e-approximate state k-design if 



E 



mm 



E 



mm 



< 



/k+d~l\ 
\ d~l ) 



(2.4) 



('^d^i ^) ^-PPsars because it is the dimension of the symmetric subspace. 

We will need a definition of an approximate unitary design and will use a slightly different form to 
the approximate state design definition above that is simpler for our purposes: 

Definition 2.6. v is an e-approximate unitary k-design if, for all balanced monomials M of degree 
< k, 



\E^M{U) -EhM{U)\ < 



d^ 



(2.5) 



Finally, we will say a few words about the efficiency of implementing an approximate unitary 
design. We would like to be able to have an e-approximate fc-design from which we can sample and 
implement the unitaries using poly(logd, fc, log 1/e) resources. Firstly, Ref. [2] provides an efficient 
construction of an e-approximate state /j-design for all k < d/2. For unitary designs, we can use the 
efficient tensor product expander construction from [13|. A {d,D,X,k) tensor product expander 
(TPE) is an ensemble of D unitaries v in dimension d with, for all p, 



U®^p{U^)®^ 



E 



H 



< ^ IIpII 



(2.6) 



where A < 1. In [13], an efficient construction is presented with D and A constant for k = 
0(log(i/loglogd). In particular, we can obtain an efficient construction for k = loiog^^ipg ^j - To 
obtain a design according to Definition 12.61 '^^ can iterate the expander: 



6 



Lemma 2.7. Iterating a {d, D, A, k)-TPE 0{klogd + \og 1/e) times gives an e- approximate unitary 
k- design. 

The slightly technical proof is in the Appendix lA.il Using the efficient TPE construction from |13j . 
we have an efficient construction of an e-approximate /c-design for k = 0(log d/ log log d). 



3 Main Technique 

The main idea in this paper can be summarised in three steps. Let / : U{d) ^ C be a polynomial 
of degree K in the matrix elements of a unitary U. Then to get a concentration bound on / when 
U is chosen from a fc-design: 

1. Find some measure concentration result for \f{U) — /x| when the unitaries are chosen uniformly 
at random from the Haar measure. Normally /i will be the expectation of /. 

2. Use this to bound the moments E|/(C/) — ^p"* for some integer m < 

3. Then use Markov's inequality and the fact that for a (approximate) fe-design the moments 
are (almost) the same as for uniform randomness. We then optimise the bound for m, which 
will often involve setting m close to the maximum, \_-^\ ■ 

We will now work through each of these steps and finish with a proof of Theorem 11.21 



3.1 Step 1: Concentration for uniform randomness 

For the first step, we will often start with Levy's Lemma. This states, roughly speaking, that slowly 
varying functions in high dimensions are approximately constant. We quantify 'slowly varying' by 
the Lipschitz constant: 

Definition 3.1. The Lipschitz constant rj (with respect to the Euclidean norm) for a function f is 

. \fiUi)-f{U2)\ 
r? = mm — r- r- -. o.l 

Then we have Levy's lemma: 

Lemma 3.2 (Levy, see e.g. [17|). Let f be an r]-Lipschitz function on U{d) with mean E/. Then 

m - m > < 4exp (-^^) (3.2) 

where Ci can be taken to be 
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3.2 Step 2: A bound on the moments 

Levy's Lemma says that / is close to its mean. This means that E|/ — E/l™ should be small. We 
will bound the moments for slightly more general concentration results: 

Lemma 3.3. Let X be any random variable with probability concentration 

F{\X - fi\>5) <Ce-''^\ (3.3) 
(Normally fi will be the expectation of X, although the bound does not assume this.) Then 



1^ / m \ m/2 

mx - < CTim/2 + l)a-'"/2 < C ( — ) (3.4) 

V2a/ 



for any m > 0. 



Proof. This proof is based on the proof of an analogous result by Bellare and Rompel [5], Lemma 
A.L 

Note that, for any random variable 1^ > 0, 

EY = F{Y > y)dy. (3.5) 
Jo 

Therefore 



j-OO 

E|x-/ir= / p(|x - /^r > x)(ix 

Jo 

POO 

= / F{\X >x^/™)dx 
Jo 

<C exp(-aa;2/™)dx 
Jo 



where in the last line we used the assumed large deviation bound Eqn. 13.31 To evaluate this integral, 
use the change of variables y = ax^l^ to get 

E|X - /^r < — a-™/2 / e-yy^l'^-^y 
2 Jo 



= Ca-'^l'^Tijnll + 1 
< 



\2a 



3.3 Step 3: A concentration bound for a fc-design 

Now we show how to obtain a measure concentration result for polynomials when the unitaries 
are selected from an approximate /c-design. We first show that the moments of \f — \x\ for / a 
polynomial are close to the Haar measure moments: 
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Lemma 3.4. Let f be a polynomial of degree K and fi be any constant. Let f = Yll=i o^i^i where 
each Mi is a monomial. Let a{f) = Then for m an integer with 2mK < k and v an 

e- approximate k-design, 



E,i/ - /.|^- < ehI/ - Mr™ + ^ (« + ImI) 



\2m 



(3.6) 



Proof. For simplicity, we assume that / and fi are real. Our proof easily generalises to the complex 
case. 

Firstly we calculate |Ei//* — Mnf^l using the multinomial theorem: 



ki+...+kt=i 



< 



E 



ki+...+kt=i 



ki,...,kt 



ki,...,kt 



... ( E^A/f 1 . . . * - Ej^Mf ^ . . . M/ 



fkl ji/rkt 



E^Mf . . . - E^Mf . . . 



~ S ^ \ki,...,kt 

ki+...+kt=i ^ 



|ai|^i...|at|'=* 



We now calculate E,^ | / — ^ 



2m. 



e.|/-mP'"-Eh|/-m 



2m| 



e.(/-m)''"-e^(/-/.) 

'£(^"){^ur-^Hf'){-t^) 



2m— i 



< 



i=0 
2m 

E 

i=0 



2m 



iE.r-EHriiM 



i 1 1 , , 1 2m— i 



siE(T)°>i^"- 

= ^(« + l^l)""■ 



□ 



Now we can simply apply Markov's inequality to prove Theorem 11.21 
Proof of Theorem Apply Markov's inequality and Lemmas 13.31 and [37 



< 



□ 
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We finish this section with two remarks. Firstly, provided a{f) (the sum of the absolute value of 
all the coefficients) is at most polynomially large in d, we can choose e to be polynomially small to 
cancel this at no change to the asymptotic efficiency. Secondly, when applying the theorem we will 
optimise the choice of m (and normally choose k = 2mK). Often a = Q{d) and the optimal choice 
of m is often Q{d) as well. However, we will not take m so large because we can only implement 
an efficient fc-design for k = O (log d/ log log c?). 



4 Application 1: Entropy of a /c-design 

We now apply the above to show that most unitaries in a /c-design generate large amounts of entropy 
across any bipartition, provided the dimensions are sufficiently far apart. This means that, for any 
initial state, for most choices of a unitary from a /c-design applied to the state, the resulting state 
will be highly entangled. We go via the purity of the reduced density matrix, since the entropy 
function is not a polynomial. 

We will call the two systems S (the 'system') and E (the 'environment') and calculate the purity 
of the reduced state. That the purity, tr (tiEUpU^)'^ , is a polynomial of degree 2 is easily seen 
by noting that the trace is linear and the reduced state is squared. However, we should check that 
there are not too many terms or terms with large coefficients. To do this, we should calculate a to 
apply Theorem II. 2 [ 

There is a general method for calculating a{f) which we will use. Write f{U) = ^•ajMj(C/) for 
monomials Mj. To evaluate a(/) = Y2i calculate f{A) where A is the matrix with all entries 
equal to 1 (so that Mi(A) = 1) and replace Qj with \ai\. Using this here we find 

a = Y^\pij\ 

\ ij 



= d'\\p\\l 
<d\ 

We now calculate the expected purity: 

Lemma 4.1. The expected purity of the reduced state is ^^^y^, where ds is the dimension of 
subsystem S and dE = d/ds is the dimension of subsystem E. 

Proof. We have 

¥.h\\^Ijs\\1 = ^h \iTJ's,sA^TEUpU^ (^iTEUpU^)\ (4.1) 
where J~SiS2 is swap acting between systems Si and 82- By linearity of the trace, we can commute 



10 



the Eh through and use [f/pC/t ® U pW] = to find 



EH\\^ps\\^ = tr 



d{d + l\ 



1 



d{d + l) 
ds + dE 



{d%ds + dEdl) 



d + 1 

Working out the higher moments in this way is difficult so we use Levy's Lemma and Lemma [37 
To use Levy's Lemma, all we have to do is find the Lipschitz constant for the purity: 

Lemma 4.2. The Lipschitz constant for purity is < 2. 



□ 



Proof. 



r] = sup ■ 



sup 



b - \\(l>s\\2\ iWi^sh + \\(Ps\\2) 



Now we use |||5'||2 — 11^1121 < ||5 — T||2 to find 

r/<sup(||V5l|2 + ||'/'5l|2) <2 

'4>,<t> 

using the fact that the purity is upper bounded by 1. □ 
Lemma 4.3. For ji = "^^^t^ and m an integer with m < k/A and v an e- approximate k-design, 



MSi^Ps) < -log2/x-a) < 



(/i(2° - 



4m 



(4.2) 



Proof. We use the fact that von Neumann entropy is lower bounded by the Renyi 2-entropy 



I.e. 



■log2llV'5|li: 



log2 



^5lli. 



(4.3) 



Then 



K{S{^s) < -log2(l + <5)^) <P.(52(V'5) < -log2(l + <5)/i) 

= K{\\i's\\2 > {i + S)fi) 

<F,(|||^5||2-/i| ><5/u) 



< 



1 



2m 



4m 



using Theorem 11.21 in the last line. 



□ 



We have written this in a more convenient form in Theorem 11.31 which is proved in the Appendix 
IA.2I This is to be compared with the fully random version Theorem 11.11 As expected, we have 
n = log2 d appearing in the exponent rather than d. Note also that our bound does not work well 
for ds ~ dE- In fact, in this case, we do not get a bound that improves with dimension. To get 
this, ds must be polynomially smaller than dE- 
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5 Application 2: /c-designs and Statistical Mechanics 



We can also apply these ideas to derandomise some of the arguments on the foundations of statistical 
mechanics in [IQ]. In this paper, the authors develop the idea that the uncertainty in statistical 
mechanics comes from entanglement rather than the traditional assumption of the principle of 
equal a priori probabilities. They consider the universe being in a pure quantum state and that 
the uncertainty in the state of a subsystem comes from the entanglement between this system and 
the rest of the universe. 

The setting is that there is an arbitrary global constraint R. Often this will be a total energy 
constraint although this is not assumed. Let the Hilbert space of states satisfying R be Tin. Then 
let the system and environment Hilbert spaces be TCs and TCe respectively. Then 

nnQHs^nE. (5.1) 

Let the dimensions be dji, ds and dE and let £r = Note that dn < dsdE, unlike in the 
above where we took d = dsdE- Normally we will have ds ^ dji. The principle of equal a priori 
probabilities says that the state of the universe is £r which implies the subsystem state is the 
canonical state, given by 

ns = tlE{£R). (5.2) 

The main result of [19] (the 'principle of apparently equal a priori probabilities') is that, for almost 
all pure states of the universe, the subsystem state is almost exactly the canonical state. 

Theorem 5.1 (Theorem 1 of [E]). For a randomly chosen state G TCr C TCs ®T~iE (^nd 
arbitrary e > 0, the distance between the reduced density matrix of the system ps = ti E{\<i)){4>\) '^''^d 
the canonical state i^s (Eqn. \ 5.2\) is given probabilistically by 



PH[\\ps-^s\\i>e + J^] <2exp(-C2dRe2) (5.3) 



where C2 = 1/(18^3) df = ^>^. 



This result gives compelling evidence to replace the principle of equal a priori probabilities with 
the principle of apparently equal a priori probabilities, but it does not address the problem of how 
the system reaches this state. It will take an extremely (exponentially) long time for the universe 
to reach a random pure state, in contrast to the observed fact that thermalisation occurs quickly. 
Here, we show that for almost all unitaries in a /c-design applied to the universe, the subsystem 
state is close to the canonical state. Since these unitaries can be implemented and sampled from 
efficiently, this means that equilibrium could be reached quickly to match observations. 

We are now ready to show that a A;-design gives a small \\ps — ^^5||i. First, we have to modify 
Lemma 13.31 slightly: 

Lemma 5.2. Let X be any non-negative random variable with probability concentration 

P(X >5 + ri)< Ce-"^' (5.4) 

where r] > 0. Then 

EX" < C ( ^ J + (2??)"^ (5.5) 
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for any m > 0. 

The proof is very similar to the proof of Lemma 13. 3[ 
Now we state and prove the main result in this section: 

Theorem 5.3. Let v he an e- approximate unitary k-design. Then 

KiWps - ns\h > s) < f (2 (-^Y' + + i(4 + 1)^-1 . (5.6) 



\ \2C2dRj \dRj '4 

In particular, with e = | (^^~") ; k < 8C2d'g, 



A\\ps-^s\\i>S)<6^-^j . (5.7) 



Again, we need ds to be polynomially smaller than dji to obtain non-trivial bounds. 



Proof. We go via the 2-norm and use Lemmas 15.21 and 13.41 
We have from Theorem 15.11 that 

rH{\\ps-^s\\i>S + v)<2e-^^''^^' (5.8) 
where r] = < Since \\ps - fish < \\Ps - ^sWi, 

P//(||p5-^^s||2 >5 + ??) <2e-^2'^«^'. (5.9) 
We now apply Lemma 15.21 to get 

EH\\ps-ns\\l'^<2^-^j +(2r/)2-. (5.10) 

So for m < k/4, using Markov's inequality and Lemma 13.41 (with /i = 0) on the polynomial 

Wps - ^sWl ■ 

KiWps - ns\U > ^) < ^ (2 {^y + (2^)'"^ + iridl + 1)'"^) ■ (5-11) 

Here, we used an estimate of a, the sum of the moduli of the coefficients: 

a < (4 + 1)2 (5.12) 

which we obtain via a similar calculation to that in Section HI 

Now we go back to the 1-norm, using \ \ps — ^s\\i < VdsWps — ^s\\2 to get 

MWps - ^sWi >S)< m\\ps - nsh > 5/^/d~s) (5.13) 
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To obtain the result in Eqn. 15.61 we just use r] < -^jj^- 

To prove the simphfied version, first use, as in Section HI that + 1)^™ < 2d^^ for m < 
This is imphed hy k < SC2<ig- We then set m = k/S to find 

P„(I|P. - > < 2 + + 2^. (5.15, 

Then, using k < 8C2(i|, with e < | (^^) , we obtain the simphfied result Eqn. 15.71 □ 



6 Application 3: Using /c-designs for Measurement Based Quan- 
tum Computing 

Here we apply our ideas to derandomise some results of Gross, Flammia and Eisert in [llj and 
Bremner, Mora and Winter in [6j. The main result in these two papers is that most states do 
not offer any advantage over classical computation when used in the measurement based quantum 
computing (MBQC) model. In MBQC, a classical computer is given access to a large quantum 
state on which it can do single qubit measurements. Some states allow for universal quantum 
computation whereas others do not add any extra power to the classical computer. These results 
are concerned with the question of characterising which states do and do not work. Showing that 
random states do not give any speed up shows that useful states for MBQC are not generic and so 
must be carefully constructed. 

While the results in these two papers are similar, we will concentrate on the methods from [llj 
since their methods are simpler to apply here. They prove their result by showing that most states 
are very entangled in the geometric measure (see Definition 16. ip . They then use this to show that 
the measurement outcomes of even the best possible measurement scheme are almost completely 
random. In fact, the state could be thrown away and the measurement outcomes replaced with 
random numbers to solve the computational problem just as efficiently. This shows that you can 
classically simulate any quantum computation that uses these highly entangled states. The measure 
of entanglement they use is the geometric measure: 

Definition 6.1. The geometric measure of entanglement of a state \^) is ]21\ [j]/ 

i?,(|M/)) = -log2SupKa|^)|2. (6.1) 

where V is the set of all product states. 

They show that any MBQC using a state |^') with £'^(1^')) = n — 0(log2n) can be efficiently 
simulated classically. They then show that (we abuse notation slightly by writing Fh for P|ii-)^,y^, 
etc.) 

Theorem 6.2 (tfl], Theorem 2). For n > 11, 

¥H{Egm)<n-2log2n-3)<e-^\ (6.2) 
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This shows that most states are useless. We derandomise this result to show that most states in 
an e-approximate (e can be taken as a constant) state n^-design have high geometric measure of 
entanglement and thus are useless in the same way. 

We could apply our technique and use Theorem [L2] but in this case, it is simpler to directly bound 
the probability using Markov's inequality. 

Lemma 6.3. 

„ 777 1 / 777 \ Tfl 

P.(K<^l^)P>'^)<(l + ^)py;;r<(l + ^)y (6-3) 

where \^) is chosen from an e-approximate state k-design u, m <k and a positive integer and \ ^) 
is any fixed state. 

Proof. We prove this bound directly using Markov's inequality: 

P^(|($|*)|2 >5)= P^(|($|^)|2'" > S"^) 



1 2m 



^ V d-l ) 

1 + e 



(1 + e)m\ ,^ . f m\m 
We now prove the main result in this section: 

Theorem 6.4. For \^) randomly drawn from an e-approximate state k-design with d = 2" 

¥^,{Eg{\^)) <n-5)< (l + e)exp2(/clog2 2A: + 4nlog2l0n- A:(5 + 4n(n-(5)). (6.4) 
In particular, for k = , 5 = 31og2 n + 5 and e = 1, 

F,y{Eg{\^)) < n-31og2n-5) < 2-n-"'. (6.5) 

We note that this bound is almost the same as in Theorem 16. 2[ It only works for slightly larger 
deviations from n, which is why we obtain a slightly better probability bound. Note also that we 
can obtain an exponential bound in n (not d = 2") because the design is exponentially large in n. 

Proof. This proof closely mirrors the proof of Theorem 2 in [llj. We use the idea of a 7-net. M-y^n 
is a 7-net on product states if 

sup _ inf 1 1 |a) — |a) 1 1 < 7/2. (6-6) 
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In it is shown that such a net exists with |A/^,n| ^ (5n/7)^"'. We then proceed by showing 
that most states in the state design have small overlap with every state in the net using the union 
bound and Lemma l6.3i Finally, since every state is close to one in the net, we can show that most 
states in the design have small overlap with every product state. 

We now formalise the above. Using Lemma 16.31 and the union bound. 



V 



Now, we need to bound 



sup Ka|^>|^>572') <|AA,,J(l + e)(^)'<(^)'\l 

,71 I \ / \ I / 



2k ^ 



F^{Eg{\^)) <n-5)=Fj -log2 sup |(a|^)P <n-S] 

= F^ I sup Ka|^)P > 2-(""^M . 
\\a)ev J 

We now claim that 

sup |(a|^')p > (5' ^ sup > (572. (6.8) 

To prove this claim, let \a) be the state that achieves the supremum on the left hand side, and let 
\a) be the state closest to it in the (5'/2-net. It is shown in [11] that this implies for any |^) 



(a|^)P - |(q|^')|2| < 572. (6.9) 



Therefore 



> 6'/2. 

This implies that the supremum over all states in the net must be at least 6' /2 to prove the claim. 
We can now finish the proof. Set 6' = 2~("~'^) in Eqn. 16.81 and use Eqn. 16.71 with 7 = 6'/2 to find 

PJ sup |(a|*)|2 > 2-("-^) ) <P^ f sup |(q|^')|2 > 2-("-'^)-i 

< (1 + e) exp2(A:log2 2k + 4nlog2 lOn - k5 + 4n{n - 5)). □ 

Combining this with the arguments of [11] shows that most states in a state n^-design on n qubits 
are useless for MBQC. This shows that even many efficiently preparable states are useless. 

7 Conclusions 

We have seen how to turn large deviation bounds for Haar-random unitaries into bounds for Re- 
designs. The main technique was applied to show that unitaries from /c-designs generate large 



16 



amounts of entanglement. Then we showed that, if the dynamics of the universe produced a k- 
design, the entanglement generated would be sufficient to reproduce the principle of equal a priori 
probabilities. Finally we showed that most states in sufficiently large state designs are useless 
for measurement based quantum computing, in the sense that computation using them can be 
efficiently simulated classically. 

However, there are other bounds for which our technique does not work. Since we cannot obtain 
exponential bounds for polynomially sized designs, our technique cannot directly derandomise some 
bounds. Some results, for example showing that the oo-norm of the reduced state of a random pure 
state is close to 1/ds [12], are proven by using an e-net of states and the union bound. Since the 
e-net is exponentially large, exponentially small bounds are required. We do not know how to apply 
our idea to results of this kind and still have k = poly(log(i). (Note that we could cope with the 
e-net in Section [6] since it was just a net on product states which is considerably smaller.) 

It is also possible that our ideas could be used to completely derandomise some constructions 
(e.g. locking [151 |9]). If we could show that unitaries drawn from a /c-design work with non- 
zero probability, and come up with an efficient sampling method, then we could obtain efficient 
randomised constructions. 
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A Appendix 



Here we present some miscellaneous proofs. 



A.l Proof of Lemma 2.7 



Proof of Lemma \2. 1\ We claim that, if for all d x d matrices 



E 



H 



2 - d^^/'^ 



\P\\2 



(A.l) 



then a is an e-approximate fc-design. To prove this claim, let m < A; and take M to be any 
balanced monomial of degree m. Write M = Up^q^ . . .Up^qm^r^s^ ■ ■ -U* ^ . Then let = 



Let Ea,k{p) 



E„ 



U^^p{U^f'\8H,k{p) 



E 



H 



and pk 



17 



T<S>k — m 

Pm ® rffc-m • Then 



> \\£cr,k{Pk) - £H,k{Pk)\\2 

i,£a-,mi,Pin) £H,m{Pm)) 
1 



(:^k—m 



\fdF- 



\£a,m(.Pm) '^_H',m(/'m) I h 



We then use the fact that the largest matrix element is upper bounded by the 2-norm. For any 
matrix A, 



\Aij\ < /j;|A,,,v|2 = ViZ4U = ||A||2. 



For us, this implies 

\{£a,m{Pm) £H,m{Pin))pi...pm,'''i...rm \ — ll^o",r)i(Pm) £H,m{Pm)\\2 

which gives 



\E^M -EhM\ < 



(A.2) 
(A.3) 



to prove the claim. 

Then we just have to show how to obtain Eqn. lA.ll from Eqn. 12.61 Iterating the TPE t times gives 



< A* 



(A.4) 



where i/* is the ensemble obtained by applying t unitaries from u. Now choose t such that A* < 

□ 



^3fe/2 • 

A.2 Proof of Theorem 1.3 

Here we prove the more convenient form of Lemma 14.31 stated as Theorem 11.31 

Proof of Theorem \1.3\ Firstly, we will write the left hand side of Eqn. 14.21 in a more useful way. 
Using ln(l + x) < x, we find 

- log2 p > logs ds - f3 
where (3 = j^^, following the notation in [16]. This means 

K{S{^s) < logs ds-a-P)< MSii^s) < - logs P - 



< 



1 ( . ( 4m \ " e , ,4 ,9™ 

+ ^id^ + P) 



(/x(2" - 1))2"^ V \CidJ 
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We now simplify the right hand side. Let 5 = 2" — 1. For ds > 2, we have /x > 1/ds- We shall also 
assume that m = k/8. This gives us (using /x < 1) 

K(S(i,s) < log. * - c - ffl < ) (4 (^) . (l 4) . (A.5) 

Now, one can easily show (e.g. by induction on n) that 

(1 + 5)" < 2 (A.6) 

for 2nS < 1. We use this for n = A;/4 and S = The condition is then k < 2d^, which we shall 

assume (we will set k = log c// log log d later). We now obtain 

MSi^s) < log2 ds-a-P)< (^^) (^A + 26^ . (A.7) 

We will now take e = 2 ( 1 , so that the two terms are the same, log 1/e is poly log d so this 
remains efficient. Now 

Assuming that 6^ > we should take k as large as possible up to ^^^r^, when the right hand 

side is maximised. We then find the result after further simplification. □ 
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